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ABSTRACT 

We present an analysis of the gamma-ray data obtained with the Large Area Telescope (LAT) onboard the 
Fermi Gamma-ray Space Telescope in the direction of SNR W49B (G43.3-0.2). A bright unresolved gamma- 
ray source detected at a significance of 38er is found to coincide with SNR W49B. The energy spectrum in 
the 0.2-200 GeV range gradually steepens toward high energies. The luminosity is estimated to be 1.5 x 
10 36 (£>/8 kpc) 2 erg s" 1 in this energy range. There is no indication that the gamma-ray emission comes 
from a pulsar. Assuming that the SNR shell is the site of gamma-ray production, the observed spectrum can be 
explained either by the decay of neutral it mesons produced through the proton-proton collisions or by electron 
bremsstrahlung. The calculated energy density of relativistic particles responsible for the LAT flux is estimated 
to be remarkably large, U e . p > 10 4 eV cm -3 , for either gamma-ray production mechanism. 
Subject headings: acceleration of particles — ISM: individual(W49B) — radiation mechanisms: non-thermal 
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1. INTRODUCTION 

Galactic cosmic rays are widely believed to be accelerated 
in supernova remnants (S NRs) through the diffusive shock 
acceleration process (e.g., Blandford & Eichler 1987). Sev- 
eral SNRs have recently been detected with the Large Area 
Telescope (LAT) onboard the Fermi Gamma-ray Space Tele- 
scope; specifically SNRs W51C , Cassiopeia A, W44, and 
IC443 dAbdo et alj|2009cl l2010allcTldl) . Except for Cas A, the 
LAT-detected SNRs are known to be interacting with molec- 
ular clouds. The GeV emission from such SNRs is expected 
to be dominated by the hadronic gamma rays due to the de- 
cay of 7T° mesons, since the ambient dense molecular cloud 
would enhance the proton-proton collisions ( Aharoni an et alj 
1994). The observed gamma-ray sources associated with 
cloud-interacting SNRs are all seen to be spatially extended in 
the LAT data. Based on the extension and its comparison with 
radio data, it is concluded that the gamma-ray emission comes 
from SNRs not from pulsars/pulsar wind nebulae (PWNe). 
The LAT spectra of these SNRs steepen above a few GeV. Al- 
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though electron bremsstrahlung cannot be ruled out, tt -decay 
emission i s the most plausib l e exp lanation for the observed 
LAT data dAbdo et alJl2009a HoTOcT) . The breaks in the ob- 
served spectra may be accounted for by an energy-dependent 
escap e of accelerated protons at SNRs (Aharonia n & Atoyanl 
1996). 

The gamma-ray measurements in the TeV range provide di- 
rect support for the acceleration of p articles up to ~ 100 TeV 
in SNR shells (Aharonian et al. 2007a[). The TeV gamma rays 
in SNR RX J1713. 7-3946, one of the most prominent exam- 
ples of TeV-emitting SNRs, can be ascribed to the decay of 
7T° m esons produced in pp collisions (e.g.. lBerezhko & Volkl 
2008) if th e average magnetic fi eld strength is larger than 
~ 15 fiG dUchivama et alj [20071) . However, the emission 
mechanism remains unsettled largely because of poorly con- 
strained physical conditions in the gamma-ray-emitting zone 
in SNR RX J1713.7-3946. Other examples are valuable for 
discriminating the origins of the gamma-ray emission. 

SNR W49B (G43. 3-0.2) has a bright radio shell and cen- 
trally peaked thermal X-ray emission. The interaction be- 
tween W49B and molecular clouds was evidenced by obser- 
vation s of mid-infrared lines from shocked molecular hydro- 
gen dReach et al.ll20p6|). H I Zeeman observ ations also sug- 
gest the interaction (Brogan & Troland 2001). Near infrared 
[Fe II] emissio n exhibits filamentar y structures, tracing ra- 
diative shocks (Ke ohane et alj|2007l) . The age of W49B is 
estimated to be in the range of ~ 1000^1000 yr dPve et alj 
1984; Hwang et al. 2000 ), and the distance is estimated to 
be 8- 11 koc dRadhakrishnan et al.lll972tlMoffett & Reynolds! 
[1991: HEogan & Trol and 2001). The radio continuum map 
shows a shell structure with a diameter of ~ 4' (~ 10 pc 
at 8 kpc). The radio flux density is 38 Jy at 1 GHz. 
The radio emissi on is linearly polarized and the spectral in- 
dex is a = 0.48 dGreenlll988l) in the frequency range 0.3- 
30 GHz, indicating a synchrotron origin. No optical emis- 
sion is detected from the source due to the severe extinc- 
tion throug h the Galactic plane. Although the ATNF pulsar 
database ( Manchest er et alj I20 05F" 4 ! lists seven pulsars with 
the spin-down luminosity > 1 x 10 34 erg s" 1 within 1 .°0 of the 
SNR position, no pulsar candidate has been reported within 
0.°4. Prior to our LAT observations, gamma-ray emission had 
not been detected in the GeV or TeV bands. 

Here we report the LAT observations in the direction of 
SNR W49B. A GeV gamma-ray source spatially coincident 
with W49B is designated as 0FGL J 1 9 1 1 .0+0905 in the initial 
source list published by the Fermi LAT collaboration, which 
includes the 205 most significant sources based on t he ob- 
servation in the first three months ( Abdo et al. 2009b). Is is 
also designated as 1FGL J1910.9+0906c in the year-1 catalog 
(1FGL catalogQ In this paper, we present a detailed anal- 
ysis of this LAT source with much longer accumulation time 
of about 17 months. This paper is organized as follows. In 
section 2, the observation and the data reduction are summa- 
rized. The analysis results for the LAT source in the direction 
of SNR W49B are reported using 17 months of the LAT data 
in section 3. In section 4, we discuss whether the gamma rays 
come from the SNR shell or a pulsar, and study the cosmic -ray 
acceleration using multi-wavelength data. 

2. OBSERVATION AND DATA REDUCTION 
s4 http://www.atnf.csiro.au/research/pulsar/psrcat 

5S available at th e Fermi Science Support Center. 
http://fermi.gsfc.nasa.gov/ssc 
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The Fermi Gamma-ray Space Telescope was launched on 
June 11 2008. The Large Area Telescope (LAT) onboard 
Fermi is composed of electron-positron pair trackers, featur- 
ing solid state silicon trackers and cesium iodide calorimeters, 
sensitive to photons in a very broad energy band (from 0.02 to 
>300 GeV). The LAT has a large effective area (~8000 cm 2 
above 1 GeV if on-axis), viewing ^2.4 sr of the full sky with a 
good angular resolution (68 % containment radius better than 
~1° above 1 GeV). The tracker of the LAT is divided into 
front and back sections. The front section (first 12 planes) has 
thin converters to improve the PSF, while the back section (4 
planes after the front section) has thicker converters to enlarge 
the effective area. The angular resolution of the back events 
is a factor of two worse than that of the front events at 1 GeV. 

The LAT data used here were collected for about 17 months 
from August 4 2008 to December 26 2009. The diffuse event 
class was chosen and photons beyond the earth zenith angle 
of 105° were excluded to minimize Earth albedo gamma rays. 

Among the standard science analysis toolfl we utilized 
gtlike for spectral fits and gtfindsrc to find a point source 
location. With gtlike, an unbinned maximum likelihood fit 
is performed on the spatial and spectral distributions of ob- 
served gamma rays to optimize spectral parameters of the in- 
put model taking into account the energy dependence of the 
point spread function (PSF). On the other hand, gtfindsrc op- 
timizes a point source location by finding the best likelihood 
for different positions around an initial guess until the con- 
vergence tolerance for a positional fit is reached. The P6_V3 
instrument response functions were used for the analyses in 
this paper. Deta ils of the LAT instr ument and data reduction 
are described in lAtwood et al.l (120091) . 

3. ANALYSIS AND RESULTS 
3.1. Detection and Source Localization 

The LAT observation revealed significant (38er) gamma-ray 
emission from the direction of SNR W49B with 17 months 
of data. Figure Q] shows LAT count maps in the vicinity of 
SNR W49B in the 2-6 GeV and 6-30 GeV bands. Only front 
events are used in the count map to achieve better angular 
resolution. The effective LAT point-spread function (PSF) is 
constructed using a spectral shape obtained through a maxi- 
mum likelihood fit (gtlike) in t he corresponding energy band 
for each count map (see jj3.31 l. The statistical and system- 
atic uncertainties in the spectral shape do not noticeably af- 
fect the PSF shape. A Spitzer near-infrared (5.8 /im) map, 
which traces ionic shocks in the SNR, is overlaid on the count 
mapQ Both count maps clearly suggest that gamma-ray 
emission comes predominantly from the SNR W49B region, 
not from a nearby star forming region, W49A. Comparisons 
between gamma-ray distributions and LAT PSFs in both en- 
ergy bands indicate that the observed gamma-ray emission 
could be consistent with a point source. 

In order to confirm the consistency with a point source, a 
radial profile of the gamma rays from the above source lo- 
cation is compared with that expected for a point source for 
front events in 2-30 GeV band as shown in Figure [2] The 
background, which is composed mainly of the Galactic dif- 
fuse emission, is subtracted. No sign of spatial extension can 
be seen in Figure [2] 

To evaluate the consistency with a point source quantita- 
tively, we compared the likelihood of the spectral fit for a 

56 The IR data are available from NASA/IPAC Infrared Science Archive. 
http://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE 



point source and an elliptical shape (3' x 4' in size; compat- 
ible with the extent of the IR image as shown in Figure [TJ 
with a uniform surface brightness. Here, we assumed a bro- 
ken pow er la w function to model the source spectrum in the 
fit (see 33.31 for details). The resulting likelihood was al- 
most the same for both cases (difference of log likelihood was 
~ 3), which means the source emission is consistent with that 
from a point source. Therefore, to simplify the analyses, the 
gamma-ray source in the SNR W49B region is analyzed as 
a point source in this paper. Assuming a point source, the 
gamma-ray source position was found to be (a, <5)=(287.°756, 
9.°096) with an error radius of 0.°024 at 95% confidence level 
using gtfindsrc, as indicated by the black circle in Figure[T] 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 




Fig. 1. — Fermi LAT count map in the vicinity of SNR W49B in 
units of counts per pixel. The pixel size is 0.°01. The LAT localiza- 
tion is represented by a black circle with a radius of 0.°024 (95 % 
confidence level) centered at (a, <5)=(287.°756, 9.°096). Cyan circles 
represent radii of the effective LAT PSF at 75%, 50% and 25% of 
the peak. Magenta and green contours indicate W49B and W49A 
in the Spitzer IRAC 5.8 /im, respectively. (Top) The count map in 
2-6 GeV is smoothed by a Gaussian kernel of a = 0.°2. (Bottom) 
The count map in 6-30 GeV is smoothed by a Gaussian kernel of 

CT = 0.°1. 



3.2. Evaluation of Galactic Diffuse Model 

Since uncertainties associated with the underlying Galac- 
tic diffuse emission are expected to be the largest system- 
atic effects for spectral analyses of the W49B source, those 
effects should be carefully evaluated. The uncertainties of 
the Galactic diffuse emission are primarily due to the im- 
perfection of the Galactic diffuse model and/or the contribu- 
tions from unresolved point sources. As a first step of the 
evaluation process, the position and energy dependences of 
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Fig. 2. — The radial profile of the LAT count map with front data 
in 2 to 30 GeV in units of counts per pixel. The pixel size is 0.°01. 
The origin of the profile is the LAT localization of the W49B source 
at (a, <5)=(287.°756, 9.° 096). The histogram shows the profile of the 
effective LAT PSF. 

the discrepancies between the observed gamma-ray distribu- 
tions and the Galactic diffuse model are studied in the re- 
gions where the Galactic diffuse emission is considered to 
be dominant around the W49B source. The normalization 
of the Galactic diffuse model is determined by running gt- 
like for a circular region with a radius of 10° centered on the 
W49B source in the energy range of 0.2-200 GeV. The posi- 
tion of the W49B source is fixed at (a, cf)=(287. 756, 9.°096) 
determined by gtfindsrc (see § 13. It . The positions and spec- 
tral shapes of all other sources are fixed at the value in the 
1FGL catalog, while the flux is al lowed to vary, except for 
PSR J 1907+06 (lAbdo et alJl2010bh. SNR W51C (lAbdo etalJ 
I2009d) and SNR W44 dAbdo et all l2010cl) which are 3, 6 
and 9 degrees away from W49B, respectively. Since these 
sources around SNR W49B are very bright as evident in 
Figure [3] we carefully evaluated spectral models for these 
sources. For this study we modeled W44 as two point sources 
at (a, <5)=(283.°89, l.°56) and (284.° 10, l.°15), to approxi- 
mately account for its angular extent, while the positions of 
the other sources are fixed at the values determined by the 
catalog. The spectral shape of these four bright sources is 
assumed to be a broken power law since likelihood tests be- 
tween a power law function and a broken power law func- 
tion favored a broken power law hypothesis at > 10ct (PSR 
J1907+06), > la (SNR W51C), > 15cr (SNR W44) confi- 
dence levels. All spectral parameters (flux, spectral break 
and spectral indices at low and high energy) are allowed to 
vary in the fitting since the spectral model is different from 
that reported in the 1FGL catalog. The Galactic diffuse emis- 
sion is modeled using "gll_iem_v02.fit". An isotropic com- 
ponent (isotropic_iem_v02.txt) is also included to account for 
instrumental and extragalactic diffuse backgrounds. The both 
background models are the standard diffuse emission models 
released by the LAT tearrQ. The normalization factors of the 
Galactic diffuse and the isotropic models are allowed to vary. 

In order to evaluate the the validity of spectra for back- 
ground models, we compared two counts spectra in the 0.2- 
10 GeV band; the spectrum expected from the models ob- 
tained by the above procedure, and the observed spectrum. 
This comparison is performed in a nearby circular region with 
a radius of 0.°5 centered on A/ ~ +2° and Ab ~ 0° from 
SNR W49B where the Galactic diffuse component is domi- 
nant. Figure [4] shows resulting fractional residuals, namely 
(observed-model)/model, as a function of energy. We fit the 
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Fig. 3.— The LAT count map above 1 GeV around SNR W49B in 
units of counts per pixel. The pixel size is 0.° 1, and Gaussian smooth- 
ing is applied with a kernel size of a = 0.°3. The W49B source, 
which is located at the center, is clearly visible. The green boxes 
(l.°0 x l.°6) represent the regions used for the evaluation of spatial 
dispersion of the difference between the Galactic diffuse model and 
the observed distribution. The cyan and white circles represent the 
regions where the flux of the Galactic diffuse model was varied to 
evaluate effects of spatial dispersion of the model. 

residuals with a cubic function as shown in Figure [4] which 
will be used to estimate the systematic error in flux d ue to 
uncertainties of Galactic diffuse model as discussed in § 13.31 
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Fig. 4. — The fractional residuals at 8 energy bins in 0.2 to 10 
GeV between the observed LAT data and the best-fit Galactic dif- 
fuse emission model in the nearby circle centered at (a, 5)=(288.°4, 
10.°2) with radius of 0P5. The fluxes of all sources included in the 
fit-model except for the Galactic diffuse component are subtracted 
from the observed data. The blue line shows a cubic function fitted 
to the residual data. 

Uncertainties of the spatial distribution of the Galactic dif- 
fuse emission are evaluated by measuring the dispersion of the 
fractional residuals in 14 regions, where the Galactic diffuse 
component is dominant (Figure[3). The regions around 4 very 
bright sources, the W49B source, PSR J 1907+06, SNR W51C 
and SNR W44, are excluded. The fractional residual for each 
region is calculated in 5 energy bands, 0.20-0.32 GeV, 0.32- 
0.50 GeV, 0.50-0.80 GeV, 0.80-1.3, 1.3-10 GeV. Figure 
shows the resulting distribution of the fractional residuals for 
14 regions in 5 energy bands. The figure shows that 68% and 
90% of the fractional residual are within 4% and 6%, respec- 
tively. To be conservative, the fractional residual of 6% will 
be used to estimate the systematic error in flux due to uncer- 
tainties of Galactic diffuse model below. 
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Fig. 5. — The histogram of the fractional residuals between the ob- 
served LAT data and the Galactic diffuse model fixed at the best-fit 
parameters determined by gtlike. The fluxes of all sources included 
in the fit-model except for the Galactic diffuse component are sub- 
tracted from the observed data. The residual was calculated in 5 en- 
ergy bands (0.20-0.32 GeV, 0.32-0.50 GeV, 0.50-0.80 GeV, 0.80- 
1.3, 1.3-10 GeV) for each region as shown in Figure[3]- 

3.3. Gamma-ray Spectrum ofW49B 

The spectral energy distribution (SED) of the source asso- 
ciated with W49B is evaluated by dividing the 0.2-200 GeV 
energy band into 1 1 energy bins, extracting data inside a cir- 
cular region with a radius of 10° centered on the W49B source 
and by using gtlike to obtain a flux value at the center of each 
bin. In each gtlike run, the W49B source, the other 1FGL 
sources, Galactic diffuse and isotropic backgrounds are fit- 
ted with their normalization free. The W49B source is fitted 
with a simple power-law function in each energy bin with its 
spectral index fixed at 2.2 below 5 GeV and 2.9 above 5 GeV 
using the fitting result in 0.2-200 GeV (see below), while the 
indices of the other sources are fixed at the values in the 1 
FGL catalog. Note that the obtained flux of the W49B source 
is insensitive to the choice of the index, if it is fixed in a rea- 
sonable range (say, 2-3). Figure [6] shows the resulting SED 
for the W49B source. 




E [eV] 

Fig. 6. — Spectral energy distribution of the W49B source measured 
with the Fermi LAT. The vertical red and black lines represent sta- 
tistical errors (la) and systematic errors, respectively. The blue line 
represents the best-fit broken power law from an unbinned likelihood 
fit in 0.2-200 GeV. 

In order to evaluate systematic effects on the SED due to un- 
certainties of the Galactic diffuse model, we varied the Galac- 
tic diffuse model used in the fit. Systematic errors due to 



uncertainties in the energy spectrum of the Galactic diffuse 
model are estimated by comparing the fit with and without the 
modification of the energy distributions of the Galactic diffuse 
model according to the curve in Figure [4] We did not modify 
the shape of the energy distribution above 10 GeV. Though 
the fractional residual intensities in Figure |4] are within our 
current understanding of the systematic uncertainties in the 
effective area, these residuals were used conservatively as un- 
certainties of the Galactic diffuse model. We obtain an esti- 
mate of uncertainties as: < 30% for below 1 GeV, < 20% 
in 1-2 GeV, and <10% above 3 GeV. Systematic errors due 
to uncertainties of spatial distribution of the Galactic diffuse 
model as shown in Figure[5]are estimated, using two modified 
Galactic diffuse models in which fluxes are varied by 6% in 
all energy bins for one of two regions with a 3 degree radius, 
a disk centered on W49B or an offset disk at (a, <5)=(285.°18, 
4.°51). The resulting systematic errors are estimated to be 
45% at 300 MeV, decreasing to 12% at 700 MeV, and <6% 
above 1 GeV for both cases. We adopt the maximum value 
among these errors at each energy bin as the systematic er- 
ror due to the Galactic diffuse model. Other systematic errors 
include uncertainties of the effective area which are 10% at 
100 MeV, decreasing to 5% at 560 MeV, and increasing to 
20% at 10 GeV and above. Total systematic errors are set by 
adding in quadrature the uncertainties due to the Galactic dif- 
fuse model and the effective area. The total systematic errors 
in each energy bin are indicated by black error bars in Figure[6] 
while statistical errors are indicated by red error bars. 

Inspection of Figure [6] suggests a spectrum steepening 
above a few GeV. We performed a likelihood-ratio test be- 
tween a power-law (the null hypothesis) and a smoothly bro- 
ken power-law functions (the alternative hypothesis) for 0.2- 
200 GeV data inside a circular region with a radius of 10° 
centered on the W49B source. The smoothly broken power- 
law function is described as 

^(^break) J ' ^ 

where photon indices Ti below the break, T2 above the break, 
a break energy Ebreak and a normalization factor K are free 
parameters. The parameter f3 is fixed at 0.05. The simple bro- 
ken power-law function is not adopted here, since the function 
cannot be differentiated at the break energy resulting in un- 
stable fit results and inaccurate error estimates. We obtained 
a test statistics of TSbpl = _ 21n(LpL/LBPL) = 22.9, which 
means a simple power law can be rejected at a significance 
of 4.4cr. The parameters obtained with the broken power law 
model are photon indices T x = 2.18 ± 0.04, T 2 = 2.9 ± 0.2, 
and £break = 4.8 ± 1 .6 GeV, with an integrated flux in 0.2-200 
GeV of (1.74 ±0.06) x 10" 7 photon cm" 2 s _1 , while the pho- 
ton index obtained with the simple power law is 2.29 ± 0.02. 
The gamma-ray luminosity in 0.2-200 GeV is calculated as 
1.5 x 10 36 (Z)/8 kpc) 2 erg s . Figure [7] shows the resulting 
fit with a broken power-law spectrum to the count spectrum 
within a radius of 0.°5 around the W49B source location. This 
underscores the importance of understanding the Galactic dif- 
fuse emission for the spectral analyses of the W49B source. 
We checked if the significance of the spectral break changes 
for different Galactic diffuse models. We found that TSbpl 
is 20.0 with the Galactic diffuse models used for evaluat- 
ing the spatial distribution uncertainties, corresponding to a 
significance of 4.1(7. TSbpl is 11.8 for the Galactic diffuse 
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model used for evaluating uncertainties of the energy spec- 
trum, corresponding to a significance of 3. Oct. Depending 
on the chosen Galactic diffuse model, the significance of the 
break ranges between 3 and 4.4ct. 




E [eV] 

Fig. 7. — The count spectrum within a radius of 0.°5 around the 
W49B source location at (a, <5)=(287.°756, 9.°096). The blue, red, 
green and black line shows best-fit model curves for the W49B 
source, the Galactic diffuse emission model, the sum of the other 
sources (all sources except for the W49B source and the Galactic 
diffuse model) and the sum of all sources included in the fit-model, 
respectively. 



3.4. Upper limit on W49A 

W49A (G43. 0+0.0) is one of the most activ e and luminous 
star fo rming regions (~ 10 7 L©) in the Galaxy (Conti & Blum 
12001) . located 0.°21 to the west of SNR W49B as shown 
in Figu re [D Its dist ance is estimated to be 1 1 .4 ± 
1 .2 kpc dGwinn et al.lfl992h . 

In this analysis, we find no gamma-ray counterpart for 
W49A. An upper limit to the GeV flux from W49A is de- 
termined by performing gtlike analysis. The model used for 
the fit includes W49A, the W49B source, all other 1FGL 
sources, Galactic diffuse and isotropic backgrounds. The 
W49A source is assumed to have uniform surface brightness 
inside a circle with radius 5'. A simple power-law function 
with its photon index fixed at 2.0 or 2.5 is used to model 
the W49A spectrum. The upper limits on the flux (0.2- 
200 GeV) obtained from the fits with the indices fixed at 2.0 
and 2.5 are 9.5 x 10~ 9 photon cm -2 s" 1 and 3.4 x 10~ 8 pho- 
ton cm -2 s" 1 at 95% confidence level, corresponding to lu- 
minosity limits of < 3 x 10 35 (D/11.4 kpc) 2 ergs" 1 and 

< 4.9 x 10 35 (D/11.4 kpc) 2 ergs" 1 , respectively. The un- 
certainties due to the Galactic diffuse model as discussed in 
33.21 have little effect on the upper limit in the case of the 
photon index 2.0, while those increase the upper limit to 

< 7.8 x 10 35 (D/11.4 kpc) 2 erg s" 1 in the case of the photon 
index 2.5. 



4. DISCUSSION 
4.1. Pulsar? 

The gamma ray emission positionally coincident with SNR 
W49B is unresolved with the LAT. This is reasonable given 
the fact that the angular extent of SNR W49B is some- 
what smaller than the effective LAT PSF The extent of 
GeV gamma-ray e mission f rom middle-aged S NRs W51C 
dAbdo et al.l2009d) and W44 dAbdo et al.l2010d) made it pos- 
sible to attribute the observed gamma-ray signals to the shells 



of these SNRs. Since this is not possible with W49B, we will 
examine a possibility that a pulsar's magnetosphere is respon- 
sible for the observed gamma-ray emission even though no 
pulsed emission has been detected with the LAT. In addition, 
no radio pulsars are found within a radius of 0.°4 around the 
LAT position of the W49B source in the ATNF catalog, while 
the LAT position is determined with 0.°024 at 95% confidence 
level. Note that we do not consider a PWN here, since the ob- 
served gamma-ray flux is very difficult to be accounted for by 
a radio-quiet PWN. 

To compare the spectral shape of the W49B sour ce with that 
of typi cal LAT pulsars in the first pulsar catalog ( Ab do et al.l 
2010e), we fit the LAT spectrum of the W49B source by a 
power law with an exponential cutoff: 



dN r 
dE F 



■—) 

^cutoff / 



(2) 



where photon index T, a cutoff energy iicutoff and a normal- 
ization factor K are free parameters. The parameters of the 
W49B source obtained by gtlike are V = 2. 10 ±0.02 and 
^cutoff = 15 ± 1 GeV. We performed a likelihood-ratio test be- 
tween a power law (the null hypothesis) and a cutoff power 
law (the alternative hypothesis) and obtained test statistics of 
TS cu toff = -21n(LpL/Lcutoff) = 27, which means that we can re- 
ject a simple power law at a significa nce of ~ 5ct. Abou t 90% 
of the 46 LAT pulsars in the catalog (lAbdo et alj|2010el) have 
r < 1.9 andZicutotf < 5.0 GeV. No pulsars exhibits ^cutoff > 6.5 
GeV among the LAT pulsars that have an error on E cuto ff less 
than 4 GeV. The LAT spectrum of W49B is different from 
what has been obtained for almost all gamma-ray pulsars so 
far. 

A pulsar may have eluded detection in X-rays due to the 
presence of bright X-ray emission from shock-heated plas- 
mas. Using 55 ks of Chandra data (PI: Holt, S. S.) we put 
an upper limit on the X-ray flux of a possible hidden pulsar 
of Fx (2-10 keV) < 6.5 x 10~ 14 erg s" 1 cm" 2 on the assump- 
tion that the pulsar spectrum is a power law with a photon 
index of 2.0. The foreground column density Nh used here 
is 6 x 10 22 cm" 2 . This corresponds to an upper limit on the 
X-ray luminosity of L x (2-10 keV) < 5 x 10 32 (D/8 kpc) 2 
erg s" 1 . The empirical correlation of the X-ray and spin-down 
luminosity of rotation powered pulsars can be written as: 



logL x = 1.341ogL sd - 15.34, 



(3) 



where Lx and L s( j are the X-ray luminosity in the 2-10 keV 
and the spin-down en ergy loss in units of erg s , respectively 
(iPossenti et alJl2002l) . This relation constrains the spin down 
luminosity of any undetected pulsars in W49B to be L s d < 1 x 
10 36 (£>/8 kpc) 2 erg s" 1 . However, the gamma-ray luminosity 
(0.2-200 GeV) of the W49B source is 1.5 x 10 36 (£>/8 kpc) 2 
erg s , which exceeds L s( j. Together with the spectral argu- 
ment, we conclude that the gamma-ray emission in the direc- 
tion of W49B is unlikely to come from a pulsar. 

4.2. Gamma Rays from the SNR Shell 

Here we consider a scenario in which the gamma-ray source 
originates in the radio-emitting shell of SNR W49B. This sce- 
nario is supported by the best-fit LAT position being coinci- 
dent with the brightest part of synchrotron radio emission as 
shown in Figure[8] The near-infrared [Fe II] emission, arising 
from warm ionized gas with a density o f order 1000 cm" 3 , co r- 
relates well with the synchrotron map (Keo hane et alj |2007). 
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We assume that the particles responsible for the LAT flux 
are distributed in a radio-emitting zone which can be char- 
acterized by a constant hydrogen density «h and magnetic 
field strength B. The volume of the emission zone is writ- 
ten as V = f(4ir/3)R 3 where / < 1 denotes a filling factor and 
R = 4.4 pc is the radius of the remnant. The radio-emitting 
material would originate in swept-up stellar wind and/or in- 
terstellar gas. We adopt the total mass contained in the zone 
as Mu = 50M Q , which would be valid within a factor of few. 
We then consider three cases: (1) nu = 10 cm" 3 and / = 0.6; 
(2) n H = 100 cm" 3 and / = 0.06; (3) n H = 1000 cm" 3 and 
/ = 0.006. Note that the constant product of friu implies the 
fixed mass in the gamma-ray-emitting region. Case (1) is con- 
sidered for a reference purpose, even though it would hardly 
explain the similarity between the synchrotron and the [Fe II] 
images. This set of parameters is more appropriate for the X- 
ray emitting gas, w hose density is estimated as n ~ 5-8 cm" 3 
(Micelietal. 2006). 




Fig. 8. — LAT source position at a 95 % confidence level (a magenta 
circle) is superposed on the Spitzer IRAC 5.8 /im ima ge. Contours 
show 20 cm radio intensity obtained from MAGPIS 1 Helfa nd et alj 

12003) . 

We adopt the follo wing form as injec tion distributions of 
protons and electrons (Abdo et al. 20093): 



Qe, P (p) = a e , P — 




-As/2 



(4) 



where po = 1 GeVc" . The indices and the break momentum 
are set to be common between el ectrons and p rotons. The 
radio synchrotron index a = 0.48 (lGreenlll988l) corresponds 
to s ~ 2. The kinetic equation for the momentum distribution 
of high-energy particles in the shell can be written as: 







dt dp 



(be, P N e , p ) + Q e , p , 



(5) 



where b ep = -dp/dt is the momentum loss rate, and Q e . p {p) 
(assumed to be time-independent) is the particle injection rate. 
To obtain the radiation spectra from the remnant, N e . p (p, To) is 
numerically calculated for To = 2000 yr. Note that energy loss 
processes such as ionization/Coulomb and synchrotron losses 
are generally not fast enough to modify the gamma-ray spec- 
trum in the LAT band. The gamma-ray emission mechanisms 
include the 7r°-decay gamma rays due to high-energy pro- 
tons, and bremsstrahlung and IC scattering processes by high- 
energy electrons. Calculations of the gamma-ray emissio n 
were done using the method described in 1 Abdo et al.l(l2009cl) . 



The large gamma-ray luminosity of L 7 ~ 1 x 10 erg s pre- 
cludes IC scattering as a do minant cont r ibuto r to the gamma- 
ray emission as discussed in lAbdo et alj d2009cl) . Specifically, 
the total energy required in electrons would be unrealistically 
large W e = J(j- l)m e c 2 N e dp ~ 1 x 10 51 erg. We shall con- 
sider the 7r°-decay and electron bremsstrahlung models to ac- 
count for the observed gamma-ray spectrum. 
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Fig. 9. — SEDs of W49B with model curves for three cases. Cases 
al, a2 and a3 represent «h = 10, 100, 1000 cirT 3 , respectively (see Ta- 
ble[TJi. The gamma-ray emi ssion is assumed to be dom inated by ir°- 
decay. The radio emission (Moffett & Reynolds 1994) is explained 
by synchrotron radiation from primary and secondary electrons. The 
dashed line in the radio band represents the synchrotron emission 
from the secondary electrons. The gamma-ray emission is modeled 
with a combination of 7r°-decay (dashed line), bremsstrahlung (dot- 
dashed line) and IC scattering (dotted line). 

The SED of SNR W49B in the radio and gamma-ray bands 
is shown in Figure [9] together with the 7r°-decay emission 
models. The radio data are modeled by the synchrotron ra- 
diation. We construct the 7r°-decay emission models for the 
different values of n H = 10, 100, 1000 cm" 3 (Table Q). Lep- 
tonic components (synchrotron, bremsstrahlung, and IC) are 
calculated assuming a e /a p = 0.01, a value similar to what is 
observed for cosmic rays at GeV energies. Note that con- 
tributions of the secondary electrons and positrons produced 
in pp collisions are small for the sets of pa rameters that we 
adopted (Tablel; see also Abdo et al. 2010c). The secondary 
synchrotron spectrum is shown in Figure [9] (a3), where its 
flux is about 10% of the total synchrotron flux at 1 GHz for 
«h = 1000 cm -3 . The contribution of the secondaries to the 
gamma-ray emission is also small, about 10% of the electron 
bremsstrahlung components for «h = 1000 cm" 3 . 

The product of «h and W p remains almost constant irrespec- 
tive of n H : n H W p ~ 10 x 10 51 erg cm" 3 . We obtain B ~ 240 [iG 
in the case of «h = 1000 cm" 3 . The SED itself can be formally 
explained in all the cases. The energy density of relativis- 
tic protons amounts to U p ~ 1 x 10 5 eV cm" 3 . This value is 
much higher than U p ~ 100 eV cm" 3 calculated f or 7r°-decay 
domin ant modeling of middle-aged SNR W51C (Ab do et alJ 
120093) . 

In FigureQj)] the gamma-ray spectrum is modeled formally 
by relativistic bremsstrahlung of electrons. The less-luminous 
7r°-decay component is also plotted using a e ja p =1. It is 
shown that «h > 100 cm" 3 is required to reproduce the ra- 
dio spectrum. If relativistic bremsstrahlung is responsible for 
the gamma rays, the ratio of the energy density of relativis- 
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TABLE 1 

Parameters of Multiwavelength models 



Parameters 



Energetics 



Model 


a e /a p 


As 


PbT 


B 


«H 


/ 


{a)W p or (b)W e 


(a)l/ p or (b)U e 


U B 








(GeVc- 1 ) 


QxG) 


(cm" 3 ) 




(10 50 erg) 


(eV cnT 3 ) 


(eV cm 4 ) 


(Case al) 7r°-decay 


0.01 


0.7 


4 


15 


10 


0.6 


11 


1.1 X 10 s 


5.6 


(Case a2) 7r°-decay 


0.01 


0.7 


4 


60 


100 


0.06 


1.1 


1.1 X 10 5 


90 


(Case a3) 7r°-decay 


0.01 


0.7 


4 


240 


1000 


0.006 


0.10 


1.0 X 10 5 


1400 


(Case bl) Bremsstrahlung 


1.0 


1.0 


4 


5 


10 


0.6 


2.6 


2.6 X 10 4 


0.62 


(Case t>2) Bremsstrahlung 


1.0 


1.0 


4 


20 


100 


0.06 


0.23 


2.3 x 10 4 


10 


(Case b3) Bremsstrahlung 


1.0 


1.0 


4 


80 


1000 


0.006 


0.016 


1.6 x 10 4 


160 



NOTE. — Seed photons for IC include infrared (kTj R = 3 X 10 eV, Ifa = 
density, U e _ p of radiating particles are calculated for p > 10 MeV c -1 . 
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Fig. 10. — Same as Fig.[9]but the gamma-ray emission is assumed 
to be dominated by electron bremsstrahlung. Cases bl, b2 and b3 
represent rin = 10, 100, 1000 cm 4 , respectively (see TableQJ. 

tic electrons to that of magnetic fields becomes very high, 
Uc/Ub > 100 (Table Q]). The energy density required (~ 
2 x 10 4 eV cm -3 ) is also much higher than U e ~ 20 eV cm -3 
calculated for bremsstrahlung dominant modeling of W5 1C. 

5. CONCLUSIONS 

We have studied gamma-ray emission in the direction of 
SNR W49B using about 17 months of data accumulated 
by the Fermi LAT. The observed energy spectrum in 0.2- 
200 GeV exhibits steepening toward high energies, although 
a simple power-law function cannot be completely ruled out 
given the uncertainties of the Galactic diffuse model. The 
luminosity is estimated to be 1.5 x 10 36 (£>/8 kpc) 2 erg s" 1 , 
which makes this source one of the most luminous gamma- 
ray sources in the Galaxy. 

The gamma-ray source is unresolved by the LAT, which is 
consistent with the angular size of SNR W49B (~ 4' in di- 
ameter) taking into account the effective LAT PSF. Assum- 
ing a point source, the source position is found to be (a, 
<5)=(287.°756, 9.°096) with an error radius of 0.°024 at 95% 
confidence level. This result clearly shows that the gamma- 
ray emission comes predominantly from the SNR W49B re- 
gion, not from a nearby star forming region, W49A. We put 



), optical (kT opl = 0.25 eV, f/ opl = 1 eV cm ), and the CMB. The total energy, W r . p and energy 

an upper limit on the gamma-ray luminosity of W49A as 
< 3 x 10 35 (D/11.4 kpc) 2 erg s" 1 at 95% confidence level. 
The gamma-ray emission in the direction of SNR W49B is 
unlikely to come from a pulsar. The gamma-ray energy dis- 
tribution is different from that observed for other pulsars with 
the LAT. In addition, no pulsed emission has been detected 
with the LAT nor are any radio pulsars known in this direc- 
tion. 

A good match between the best-fit LAT position and the 
brightest part of synchrotron radio emission suggests that the 
gamma-ray source originates in the shell of SNR W49B. The 
LAT spectrum can be formally explained either by 7r°-decay 
gamma rays or by electron bremsstrahlung. For both cases, 
the calculated energy density of relativistic particles is evalu- 
ated to be very high, U e . p > 10 4 eV cm -3 . 
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